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A recently developed renormalization approach is used to study the electron-phonon coupling in 
many-electron systems. By starting from an Hamiltonian which includes a small gauge symmetry 
breaking field, we directly derive a BCS-like equation for the energy gap from the renormalization 
approach. The effective electron-electron interaction for Cooper pairs does not contain any singular- 
ities. Furthermore, it is found that phonon-induced particle-hole excitations only contribute to the 
attractive electron-electron interaction if their energy difference is smaller than the phonon energy. 
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I. INTRODUCTION 

The famous BCS-theory [1] of superconductivity is essentially based on the analysis of attractive interactions 
between electrons of many-particle systems [2] . As was pointed out by Frohlich [3] such an interaction can result from 
an effective coupling between electrons mediated via phonons. The recent discovery of superconductivity in magnesium 
diboride MgB2 [4] below a rather high T c of about 39 K has attracted again a lot of interest on this classical scenario 
of a phonon-mediated superconductivity However, the electron-electron interaction derived by Frohlich [3] contains 
some problems. There are certain regions in momentum space where the attractive interaction becomes singular and 
changes its sign due to a vanishing energy denominator. 

Recently, effective phonon-induced electron-electron interactions were also derived [5,6] by use of Wegner's flow 
equation method [7] and by a similarity renormalization proposed by Glatzek and Wilson [8,9]. The main idea of 
these approaches is to perform a continuous unitary transformation which leads to an expression for an effective 
electron-electron interaction which is less singular than Frohlich's result [3] . 

Recently, we have developed a renormalization approach which is based on perturbation theory [10]. This approach 
resembles Wegner's flow equation method [7] and the similarity renormalization [8,9] in some aspects. Therefore, the 
investigation of an effective phonon-induced electron-electron interaction is very useful to compare the three methods 
in more details. Therefore, in this paper we directly diagonalize the classical problem of interacting electrons and 
phonons by use of the new renormalization technique [10]. the Hamiltonian is given by 

H = £k C L C k,<x + U * b l b 1 + H l 9 * C L C (k+q),^q + 0* c{k+q) l( r C k>q) , (1) 

k,cr q k,q,CT 

which will be used to describe superconducting properties. In (1) c k a and Ck,<r are the usual creation and annihilation 
operators for electrons with wave vector k and spin a. & q and 6 q denote phonon operators with phonon energies w q . 
The electron excitation energies £k are measured from the chemical potential /i. 

The paper is organized as follows. In the next section we briefly repeat our recently developed renormalization 
approach [10]. In Sec. Ill this approach will be applied to the electron-phonon system (1) in order to derive a BCS-like 
gap equation. Furthermore, effective electron-electron interaction derived in this framework will be compared with 
the results from former approaches [3,5,6]. Finally, our conclusions are presented in Sec. IV. 

II. PROJECTOR-BASED RENORMALIZATION METHOD (PRM) 

The PRM [10] starts from the decomposition of a given many-particle Hamiltonian H into an unperturbed part Ho 
and into a perturbation Hi 



H = H + eHi =: H(e). (2) 

We assume that the eigenvalue problem Ho\n) = E^\n) of the unperturbed part Hq is known. The parameter e 
accounts for the order of perturbation processes discussed below. Let us define projection operators P> and Qa by 
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Pa and Qa are super-operators acting on usual operators A of the unitary space. Here, Pa projects on that part of A 
which is formed by all transition operators |n)(m| with energy differences \En^ — -Em | less or equal to a given cutoff 
A. The cutoff A is smaller than the cutoff A of the original Hamiltonian H. Qa is orthogonal to Pa and projects on 
high energy transitions larger than A. 

The aim is to transform the initial Hamiltonian H into an effective Hamiltonian H\ which has no matrix elements 
between eigenstates of H.q with energy differences larger than A. H\ will be constructed by use of an unitary 
transformation 



H x = e Xx He- x > 



(5) 



Due to construction the effective Hamiltonian H\ will therefore have the same eigenspectrum as the original Hamil- 
tonian Tt. The generator X\ of the transformation is anti-Hermitian, X\ = —X\. To find an appropriate generator 
X\ we use the condition that H\ has no matrix elements with transition energies larger than A, i.e., 
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has to be fulfilled. By assuming that X\ can be written as a power series in the perturbation parameter e 



X X = ex*? 



■fx® 



the effective Hamiltonian H\ can be expanded in a power series in e as well 
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Now, the high-energy parts parts Q\X^ of X\^ l> can successively be determined from Eq. (6) whereas the low-energy 



can still be chosen arbitrarily. In the following we use for convenience Pa X 
for the effective Hamiltonian H\ up to second order in Hi follows 
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so that 



H\ = Ho + P\Hi 
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Here e was set equal to 1. The quantity Lo in (9) denotes the Liouville operator of the unperturbed Hamiltonian. 
It is defined by Lo = [Ho, -A] for any operator A. Note that the perturbation expansion (9) can easily be extended 
to higher orders in e. One should also note that the correct size dependence of the Hamiltonian is automatically 
guaranteed due to the commutators appearing in (9). 

Next, let us use this perturbation theory to establish a renormalization approach by successively reducing the 
cutoff energy A. In particular, instead of eliminating high-energy excitations in one step a sequence of stepwise 
transformations is used. Thereby, we obtain an effective model which becomes diagonal in the limit A — > 0. In an 
infinitesimal formulation, the method yields renormalization equations as function of the cutoff A. To find these 
equations we start from the renormalized Hamiltonian 



H\ = Ho.x + Hi. 



(10) 



after all excitations with energy differences larger than A have been eliminated. Now we perform an additional 
renormalization of Hx by eliminating all excitations inside an energy shell between A and a smaller energy cutoff 
(A — AA) where AA > 0. The new Hamiltonian is found by use of (9) 



H(A-AA) — P(a-aa)H\ — 2"^( A ~ AA ) 
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Here, Lo,a denotes the Liouvillc operator with respect to the unperturbed part Wo. a of the A dependent Hamiltonian 
Ha- Note that the flow equations derived from Eq. (11) will lead to an approximative renormalization of Ha because 
only contributions up to second order in Hi, a are included in Eq. (11). For a concrete evaluation of Eq. (11) it is 
useful to divide the second order term on the r.h.s into two parts: The first one connects eigenstates of Ho, a with the 
same energy. This part commutes with Ho, a and can therefore be considered as renormalization of the unperturbed 
Hamiltonian 
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In contrast the second part connects eigenstates of Ho, a with different energies and represents a renormalization of 
the interaction part of the Hamiltonian 



Hi, (A-AA) - P(A-AA)Hi,A — — P(A-AA) (P(A-AA)Hi,a) 
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Note that for small AA only the mixed term, i.e., the first part on the r.h.s. of Eq. (13), contributes to the renormal- 
ization of Hi \ 



Hi, (A-AA) — P(A-AA)Hi,A ~ — P(A-AA) 
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In the limit AA — > 0, i.e. for vanishing shell width, equations (12) and (14) lead to differential equations for the 
Hamiltonian as function of the cutoff energy A. The resulting equations for the parameters of the Hamiltonian are 
called flow equations. Their solution depend on the initial values of the parameters of the Hamiltonian. Note that 
for A — > the resulting Hamiltonian only consists of the unperturbed part W(a^o) so that an effectively diagonal 
Hamiltonian is obtained. 



III. APPLICATION TO THE ELECTRON-PHONON SYSTEM 



In this section we apply the renormalization approach discussed above to the system (1) of interacting electrons and 
phonons. The aim is to decouple the electron and the phonon subsystems. The Hamiltonian (1) is gauge invariant. In 
contrast, a BCS-likc Hamiltonian breaks this symmetry. Thus, in order to describe the superconducting state of the 
system, the renormalizcd Hamiltonian should contain a symmetry breaking field. Therefore, our starting Hamiltonian 
Ha reads 

n x = n , x + n hX , (15) 

after all excitations with energies larger than A have been eliminated, where 

H ,A = Y £k C k,<x C k,<x + Y W 1 b i b 1 ~ Y ( Ak > A C k,T C -k,l + A k,A C-k,l c k,T) + Ca, (16) 
k,(j q k 

Hi, A = PaHi = P A Y (0qCk lCT C (k+q), CT &q + •< C (k + q) l( x C k,^q) ■ (17) 
k,q,<T 

The 'fields' A^a and A£ A in Ho, a couple to the operators |C^ k ^ and c_k,j.Ck,t and break the gauge invariance. 
They will take over the role of the superconducting gap function but still depend on A. The initial values for Ak.A 
and the energy shift C\ are those of the original model 

Ak,(A=A) = 0, C (A= A) = 0. (18) 

Note that renormalization contributions to the electron energies £k, the phonon energies w q , and the electron-phonon 
interactions g q have been neglected in (15). Also, additional interactions which would appear due to renormalization 
processes have been omitted. Let us first solve the eigenvalue problem of Ho, a- For this purpose, we perform a 
Bogoliubov transformation [11] and introduce new A dependent fermionic quasi-particles 
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a l.x = K,\ c l,i - v L\ c -^i, ( 19 ) 



k,A = U k,A c -k,J. + w k : ACk,T 



where 




(20) 



7io,A can be rewritten as 

W ,a = E KX (4 iA «k, A + ^, a A,a) + 1] (£k - ^k,A) + w q 6 q fe q + c a (21) 



where the fermionic excitation energies are given by \ = y e^. + | Ak.A| 2 - 

Next let us eliminate all excitations within the energy shell between A and (A — AA) by applying the renormalization 
scheme of section II. We are primarily interested in the renormalization of the gap function Ak,A- Note that for this 
case we have to consider the renormalization contribution given by (14) 

(A, AA) := W li(A _ AA) - P ( a-aa)Wi,a- (22) 

The reason is that the renormalization (12) of Ho, a only gives contributions which connects eigenstates of with 
the same energy. Thus, this part only changes the quasiparticle energies from £k.A to -EwA-AA)- In contrast, the 
renormalization (22) changes the relative weight of the operator terms in (16) and is exactly the renormalization 
needed to describe the flow of Ak,A which will be discussed in the following. 

There is no principle problem to evaluate the renormalization contributions (14). First, one expresses the creation 
and annihilation operators by the quasiparticle operators (19) and uses the relation Lo^cut x = E-^ xa^ , and an 

equivalent relation for /?£ \ to evaluate the denominator in Eq. (14). Then the quasi-particle operators (18) have 
to be transformed back to the original electron operators. The main reason for this procedure is the fact that 
the Bogoliubov transformation (19), (20) depends on the cutoff A. Thereby, a lot of terms arise which contribute 
to the renormalization (14). For convenience, we evaluate the denominator in Eq. (14) by use of the assumption 
e k 3> |Ak,A| 2 - As it turns out, the resulting flow equations still contain sums over k. Note that the approximation 
used is valid for most of the k dependent terms (except for those k values close to the Fermi momentum). Thus, 
it does not strongly affected the renormalization contributions. The two operator expressions contributing to the 
commutator in (14) are given in this approximation by 



P(A-AA)' 



)ftx,A= Yl 6 t( A ~ AA ) ~ | £ k-£(k+q) +W q |] {gqC k ,<T C (k+q),^q + h.C.j, (23) 
k,q,(T 

j~ Q(a— aa)^i,a = E / 9k ; q(A,A + A i {^ c l^A - h. c .} (24) 

L 0,A £k - £(k+q) + Wq I- 



where 

(56 k ,q(A, AA) = 6 [|e k - e {k +q) + w q | - (A - AA)] - 6 [|e k - £(k+ q ) + w q | - A] (25) 

describes the restriction to excitations on the energy shell AA. We are not interested in the renormalization of the 
phonon modes. Therefore all contributions including phonon operators are neglected. By using (23) and (24) we then 
find from (14) 
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AwaA , p |g q | 2 £6 k , q (A, AA)9 [(A - AA) - \ey - £(k , +q) 
d/ti(A, AA) = -F(a_aa) >, ; 
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k,k',q,o",o"' 

X 



{ C (k+q),<x Ck ' ffC k'.<x' C (k'+q),<x' + h.cj . (26) 

Iii the following we restrict ourselves to renormalization contributions which lead to the formation of Cooper pairs. 
Consequently, the conditions k' = — (k + q) and a' = —a have to be fulfilled so that 

- hm = (27) 

aa^o AA v ; 

= 9 [ A ~ 2 l £k ~ e ( k +q) t A - l e (k+q) - £ k + w q |] ^ (|ek - £(k+ q ) + w q | - A) x 

k,q,fT 

M 2 



£k — £ (k+q) 



{ C (k+q)^ ct -(k+q)^<T C -k,- CT Ck, CT + h.c} 



results from (26). Here, 6 [A — 2 |e k — £7 k+q ) ] is due to the projector operator P(a-aa) m (26). Note that the 
differential expression on the l.h.s. of (27) is different from the differential dHi t \/dX. This follows from the definition 
of 5Hi(\, AA) = Wi,(A-AA) — P(a-aa)^i,a which differs from AWi(A, AA) = Wwa-aa) — ^i,a due to the second 
term. Next we can simplify the 9-functions in Eq. (27) by discussing e k > £( k + q ) and £( k + q ) > £ k separately. There 
are no contributions from the latter case. For e k > £( k + q ) the contribution from the first term in the curly bracket in 
(27) and from its conjugate can be combined. By exploiting the 0-functions the result can be rewritten as 

(Wi(A,AA) sT' Xn I, >\ Iffql 2 © [w q - |e k -e (k + q) |] 
- hm = > s (leu- e (k+Q) \ + lo„ - X) i — ! — i — riLLL x 

AA^O AA ^ Vl ( 1 | £k - £ (k + q)|+^q 

XC (k+q),o- C -(k+q),-o- C -k--cr c k,cr- (28) 

Here, we have assumed g q = .9- q . Eq. (28) describes the renormalization of the A dependent Hamiltonian TL\ with 
respect to the cutoff A. Next we use (28) to derive flow equations for the parameters A kj A and C\. For this purpose, 
a factorization with respect to the full Hamiltonian TL is carried out. The final flow equations read 



dA ki A v^A/l I , ^ Iffql 2 [ w q ~ | £ k - £( k+q )|] / > 

—77— - - 2 Z^ d \ £k _£ (k+q) +^ q - Aj r r— (c_ (k+qU C (k+q) , T ) , 

aA a I k _ £ (k+q)| + w q 



(29) 



7f-E<4A,>^- « 

k 

Note that in contrast to (28), Eqs. (29) and (30) are differential equations with normal derivatives of A k> A and C\. 
This fact can be explained as follows: As discussed above, the difference between the expressions 8TCi(X, AA) = 
Wi,(a-aa) - P(a-aa)'Wi : a and A7Yi(A, AA) = Wi,(a-aa) - ^i,a is given by the quantity Q(a-aa)Wi,a which consists 
of all matrix elements of TLi.x between eigenstates of TLo.x with energy differences between (A — AA) and A. We 
arc interested in the new Hamiltonian H(a-aa) = P(a-aa)^(a-aa) which only contains transition operators between 
states with energy differences smaller than (A — AA). Therefore, all renormalization contributions which lead to matrix 
elements with energy differences larger than (A — AA) arc not relevant. Thus, we obtain differential equations for the 
parameters of W(a-aa)- 

Note that the factor 2 in front of (29) is due to the sum over a in (28). The expectation values (...) in (29) and 
(30) are formed with the full Hamiltonian TL and are independent of A. The flow equations can be easily integrated 
between the lower cutoff (A — > 0) and the cutoff A of the original model. The result is 

A A I oV" l gq | 29 l £k ^ £ (k+g)|] / „ \ fou 

A k = A k , A + 2 ^ r : r— ( C -(k+q),i C (k+q),T/ . ( 31 ) 

q |£k - £(k+q)| + 

C=C A + ^(c kiT cL ka )(A k -A k!A ) (32) 

k 
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where a short hand notation for the desired values of Ak,A and Ck.A at A = was introduced: Ak = Awa^o)) 
C = C(a_»o)- Note that Ak and C only depend on the parameters of the original system (1) and on A^a and C\. 
The initial conditions (18) for Ak,A and Ca will be used later. For A — > the renormalizcd Hamiltonian Ti = 7i(A_»o) 
reads 



Ti 



e k c k ,<r c k^ + w i b l b <i ~ E (^ k c M c -k,i + A k c- k a c k,t) + c . 



(33) 



k.tr q k 

Ti can easily be diagonalizcd by a Bogoliubov transformation according to (19) and (20) 

n = J2^ (4"k + PJ*) + ( £k - ^ k ) + E ^ b l b i + 6 ( 34 ) 

k k q 

where E k = £ k ,(A^o), a k = a k) (A— o)> and /3 k = /?k,(A^o)- 

Finally, we have to determine the expectation values in (31) and (32). Since Ti emerged from the original model Ti 
by an unitary transformation, the free energy can be calculated either from Ti or from Ti 
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where (34) was used. The required expectation values are found by functional derivative 
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Here, /(-E/ k ) denotes the Fermi function with respect to the energy E k . If we neglect higher order corrections, Eqs. (31) 
and (32) can be rewritten as 
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where the initial conditions (18) were used. Note that Eq. (37) has the form of the usual BCS-gap equation. Thus, the 
term inside the brackets {. . . } can be interpreted as the absolute value of the effective phonon-induced electron-electron 
interaction 



k,q 



2| ffq | 2 e[ W q- 


£k - £(k+q) 


] 




£ k - £(k+q) 


+ W q 



(39) 



for the formation of Cooper pairs. In contrast to the usual BCS-thcory, in the present formalism both the attractive 
electron-electron interaction (39) and the gap equation (37) were derived in one step by applying the rcnormalization 
procedure to the electron- phonon system (1). 
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Let us now compare the induced electron-elctron interaction (39) with Frohlich's result [3] 



K" = 7 -■ (40) 

Note that (40) contains a divergency at |ek — £(k+q)| — w q . Furthermore, this interaction becomes repulsive for 
|-k — -(k+q)| > w q- Thus, a cutoff function O [w q — |ek — £(k+q)|] for the electron-electron interaction (40) is intro- 
duced by hand in the usual BCS-theory to suppress repulsive contributions to this interaction. In contrast our result 
(39) has no divergency and is always attractive. Furthermore, the cutoff function in Eq. (39) shows that the attractive 
interaction results from particle-hole excitations with energies |£k — £(k+ q )| < w q . This result directly follows from 
the renormalization process. 

Recently, Mielke [6] obtained a A-dependent phonon-induced electron-electron interaction 

V^XX = - | 2|gq ' 2 , e(|e k+q - e k | + cj q - A) (41) 

|£k - £ (k+q) I + W q 

where in (41) for convenience the A-dependence of the electron and the phonon energies is suppressed. Apart from 
the A-dependent 9-function, the main difference to our result (39) is that the cutoff function [w q — |ek — e (k+q)|] 
is not present in (41). This difference may result from Mielke's way of performing the similarity transformation [8,9] 
of the electron-phonon system (1). The similarity transformation is based on the introduction of continuous unitary 
transformations and is formulated in terms of differential equations for the parameters of the Hamiltonian. Like in 
our approach (see Sec. II) also the similarity transformation leads to a band-diagonal structure of the normalized 
Hamiltonian with respect to the eigen representation of the unperturbed Hamiltonian. Due to the renormalization 
processes also new couplings occur. Mielke has first evaluated the phonon-induced electron-electron interaction (41) 
by eliminating excitations with energies larger than A. In particular, for an Einstein model with dispersion-less 
phonons of frequency loq the interaction becomes independent of A if A is chosen less than too- Of course, for this case 
and assuming £k+ q — £k| < w q Mielke's result (41) and the result (39) become the same. Keeping this interaction 
and neglecting at the same time the remaining part of the electron-phonon interaction with energies smaller than 
A a BCS-like gap equation was derived by Mielke [6]. The main difference to our result (39) is the absence of the 
cutoff- function [w q — |ek — £(k+q) |] i which demonstrates that only particle-hole excitations with energies less than 
ujq participate in the attractive interactions. However, note that by setting A = a finite value of the interaction 
remains in (41). This interaction is non-diagonal in the unperturbed Hamiltonian as used by Mielke, H^ 1 ^" = 

J2k a e k,AC k g-Ck^ + E q w q^q- This seems to be a contradiction to the allowed properties for operators at A = 
which should commute with the unperturbed Hamiltonian. 

Finally, by use of Wegner's flow equation method [7], Lcnz and Wegncr obtained the following phonon-induced 
electron-electron interaction 

T/ Lcnz/Wcgncr ^Jggj LOg 

K k,q - r -I 2 , 

[£k-£(k+q)J 

which is attractive for all k and q [5]. The result (42) is similar to (40) if w q > £k — £(k+q)| is fulfilled. In contrast 
to our result (39), the interaction (42) remains finite even for |e k — £(k+q)| > ^q- Wegner's flow equation method 
[7] as well as the similarity transformation [8,9] are based on the introduction of continuous unitary transformations. 
Both methods arc formulated in terms of differential equations for the parameters of the Hamiltonian. However, they 
differ in the generator of the continuous unitary transformation. This leads to the different results (41) and (42). A 
detailed comparison of the flow equation method and the similarity transformation can be found in Ref . [6] . 



IV. CONCLUSION 



In this paper we have applied a recently developed renormalization approach [10] to the 'classical' problem of 
interacting electrons and phonons. By adding a small field to the Hamiltonian, which break the gauge symmetry, we 
directly derive a BCS-like gap equation for the coupled electron-phonon system. In particular, it is shown that the 
derived gap function results directly from the renormalization process. The effective phonon-induced electron-electron 
interaction is deduced from the gap equation. In contrast to the Frohlich interaction [3] no singularities appear in 
the effective interaction. Furthermore, the cutoff function which is included in Frohlich's result by hand to avoid 
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repulsive contributions to the electron-electron interaction follows directly from the renormalization procedure. This 
means that phonon-induced particle-hole excitations only contribute to the attractive electron-electron interaction if 
their energies are smaller than the energy of the exchanged phonon. 
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